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Abstract 

Magnetoelectric coupling terms are derived in piezoelectric/magnetostrictive (multiferroic) thin 
film heterostructures using Landau-Ginzburg free energy expansions in terms of three order pa- 
rameters: strain, magnetization and electric polarization. Strain is eliminated using a particular 
set of interface boundary conditions. Then, a general effective medium method is used to calculate 
the ferromagnetic resonance frequency in a BaTiOs/NiFe2 04 superlattice. This method differs 
from existing methods for treating magnetoelectric heterostructures since the magnetic and elec- 
tric dipolar fields are not assumed constant but vary from one film to another. The ferromagnetic 
resonance frequency shift is calculated as a function of applied electric field and is compared to 
some experimental results. 
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I. INTRODUCTION 



Today there is much interest in magnetoelectric heterostructures due to the possibility 
of creating multifunctional devices where magnetic properties can be controlled with an 
electric field or vice- verse For example, an applied electric field has been shown to 
shift ferromagnetic resonance (FMR) frequencies {2]. This has potential application in fast 
microwave signal processing devices with lower energy consumption than current devices. 
The most promising magnetoelectric heterostructures are those containing magnetostrictive 
and piezoelectric phases because an interface strain-mediated magnetoelectric (ME) coupling 
can occur that is orders of magnitude bigger than in single-phase multiferroic materials (for 
a review, see Ref. [lj]). 

Recently we presented an effective medium method (EMM) for calculating susceptibilities 
and resonant frequencies in a heterostructure comprising alternating multiferroic (ferroelec- 
tric and antiferromagnet) and ferromagnetic thin films |3|. This EMM is based on previous 
work on dielectric (4) and magnetic {5, heterostructures respectively in the long wave- 
length electrostatic/magnetostatic limit. Only dipolar coupling was considered between the 
films. Since strain-mediated magnetoelectric coupling in heterostructures appears to be very 
important, in this paper we present a theory to calculate the strain- mediated coupling so 
that the EMM can be applied and FMR frequency shifts can be predicted as a function of 
applied electric field. 

Existing methods to calculate the high frequency magnetic and magnetoelectric suscep- 
tibility of thin film magnetoelectric composites have focused on using a simplified EMM 
and solving magnetization equations of motion with a phenomenological magnetoelectric 
coupling term |2j. The coupling term takes the form of an energy density 

&me = B ikn EiM k M n + bi jkn EiEjM k M n , (1) 

where E is the applied electric field that is assumed constant throughout the heterostructure 
and M is the magnetization in the magnetostrictive phase. The magnetoelectric coupling 
constants B ikn and h^ kn can be expressed as a combination of elastic, piezoelectric and 
magnetostrictive constants of the individual phases and are derived by considering the strain 
boundary conditions at the interfaces and by considering the symmetries of the individual 
phases [7|. Bichurin et al. have also solved simultaneous magnetization and elastic equations 
of motion and found enhanced magnetoelectric effects near mechanical resonances [8]. 



Our EMM differs from that of Bichurin et al. in several ways. Firstly, both magnetization 
and dielectric equations of motion can be written and so the magnetic, magnetoelectric and 
electric susceptibility may all be solved for, as shown in Ref. [3J. Secondly, our effective 
medium method takes into account the interface boundary conditions for the dipolar mag- 
netic and electric fields and these fields are not assumed constant through the heterostruc- 
ture, as done in other works This means that demagnetizing and depolarizing 

effects are included in our method and different resonant frequencies will be calculated as 
compared with simplified EMMs, as discussed in Ref. [3]. 

Interface strain-mediated magnetoelectric terms like Eq. (TTJ) can be added to the mag- 
netic and dielectric energy densities of the heterostructure and dipolar boundary conditions 
together with equations of motion can be solved to find the susceptibility tensor. However, 
by writing the magnetoelectric energy in terms of electric field E rather than electric polar- 
ization P, electric fields can drive magnetization but magnetic fields cannot drive electric 
polarization. We therefore aim first in Section [XT] to write a strain- mediated magnetoelectric 
energy density in terms of M and P. Writing an energy density is not strictly valid since 
magnetization and electric polarization are defined in different phases of the heterostructure. 
But the EMM works by averaging the properties of the films so both the magnetization and 
polarization may be thought of as belonging to a single effective material, therefore this 
formalism does not present a problem. 

Previously, studies have used one of two methods to derive magnetoelectric coupling in 
heterostructures. Either linear constituent equations are considered that relate strains to 
magnetic and electric fields 0, In], U| or Landau free energy expansions are made in terms 



of magnetization, polarization and strain order parameters 12Nl4|. In Section [IT] the latter 
method will be used to derive the magnetoelectric coupling. The full free energy expansions 
for both films will be made, which leads to very complicated equations for the magnetoelectric 
coupling once strain is eliminated from the equations. This is in contrast to previous works 
that have considered the Landau expansion for only the piezoelectric film, with a misfit 



strain that comes from the attached magnetostrictive film |l2l. [l4]. or that have considered 
;he expansion for only the magnetostrictive film, with a misfit strain from the piezoelectric 
13] ]. The calculation in Section HT1 therefore represents an extension of these works. 

In Section II III it is then shown how the derived strain- mediated magnetoelectric coupling 
alters the FMR frequency and also may dynamically couple magnetic and electric excita- 
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tions. This is done by calculating just the zz-component of the high-frequency magnetic 
susceptibility using the EMM presented in Ref. [3]. The FMR frequency shifts are compared 
to some existing experimental results. 

It is hoped that the general methods presented in this work will enable experimentalist to 
estimate the resonant frequencies or the strain-mediated magnetoelectric coupling strength 
in piezoelectric/magnetostrictive heterostructures. 

II. STRAIN-MEDIATED MAGNETOELECTRIC COUPLING FOR A CUBIC 
SYSTEM 

The calculation for the magnetoelectric coupling in the piezoelectric/magnetostrictive 
heterostructure will proceed as follows. Landau-Ginzburg expansions for the free energy are 
written in terms of the order parameters: electric polarization Pj, magnetization Mj, and 
strains in the respective phases, uf) and (where i,j = x,y } z). Then strain boundary 
conditions are imposed between the phases, the free energy is minimized with respect to 
strain, and finally an effective magnetoelectric coupling is derived for the heterostructure. 
Note that all the order parameters are considered constant within a phase. This is a good 
approximation for sufficiently thin films. So-called "entire-cell" effective medium methods 



15l . 1 1 €31 ] are better adapted to treating thicker films where this assumption breaks down. 

Both the magnetostrictive (denoted with a superscript m) and piezoelectric (jp) free en- 
ergies are written assuming that a cubic symmetry exists. While cubic symmetry often 
exists in the paramagnetic or paraelectric phases of these materials, usually the transition 
to magnetic/electric order is accompanied by a lowering of the crystal symmetry. Therefore, 
this treatment may at first seem like a crude approximation below the Curie temperature. 
In fact, this is not the case and the lowering of the crystal symmetry falls out naturally in 
the subsequent calculation of the strains. For example, it can be shown that assuming the 
electric polarization P is along the (111) direction in a bulk piezoelectric material, then the 
strains minimizing the free energy correspond to an elongation of the cubic cell along (111), 
in other words, a rhombohedral distortion or a transformation to space group R3c. (See 



Ref. [17[ for a more detailed discussion.) The situation is more complex for thin films that 
are mechanically coupled to other materials; very different distortions may correspond to 
the ground state as compared with bulk crystals (see, for example, Ref. [181]). 



4 



A. Landau free energy expansions 



We start by quoting the parts of the free energy expansions for the magnetostrictive and 
piezoelectric phases that do not depend on strain. The magnetic part of the free energy 
density for the cubic magnetostrictive phase is given up to fourth order, for simplicity, by 

Q: 

^fm(Oi) = K(a 2 x a 2 y + a 2 y a 2 z + a 2 z a 2 x ) - fi M H ■ a, (2) 

where K is the anisotropy constant, Mq is the saturation magnetization, H is the applied 
magnetic field and a = M/M is the magnetization vector normalized to unity or the 
magnetization direction cosines. The magnitude of the magnetization does not vary (at 
constant temperature) and therefore terms such as a 2 and af are ignored here. The exchange 
energy is ignored as the magnetization is assumed to be all aligned within a film. Any 
dynamic motion (considered in Section IHIj) involves the magnetization in the whole film 
moving coherently. 

In Eq. (J2J) the demagnetizing free energy is not included. For a thin film it is equal 
to ^iiqMqO, 2 but for a heterostructure it is typically less than this value and takes a non- 
analytic form [3]. Here it is noted that, for the material parameters used in Section Til El the 
demagnetizing energy does not alter the equilibrium magnetization direction, as discussed 
later. The effective medium method, that will be used in Section II III to find the resonant 
frequencies, calculates the dipolar energy in an implicit and elegant way, without recourse 
to a computationally intensive method such as a dipole sums method, and obtains very 
accurate results 16|. Therefore the demagnetizing effects will be considered properly in 
Section IIII1 

Most piezoelectric materials (for example, BaTiC>3, PZT and BiFeOs) undergo a first- 
order phase transition and therefore the free energy must have terms at least up to sixth 
order in Pj. Also, the magnitude of the electric polarization can vary strongly as a function 
of applied field E. The electric part of the free energy density for the cubic piezoelectric 
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phase is [20|: 

F FE (P 4 ) = ai {P 2 x + P 2 + P 2 Z ) + a n (P* + P y 4 + Pt) 
+ a l2 {P 2 x P 2 + P 2 P 2 + P 2 P 2 x ) 
+a m (J* + Pi + Pf) + u l2 ,P 2 x P 2 P 2 

+"112 [P 2 x (Py 4 + P?) + Py{Pt + P') 

+P 2 AP' + Py)] ~E P, (3) 

where the as are dielectric stiffness coefficients, measured at constant strain. A recent paper 
suggests that terms up to Pf are required in order to accurately model phase transitions in 
BaTi03 [21]. However, in the present work the temperature dependence of the polarization 
is not being studied and in fact terms of higher order than P 2 have little affect on the 
magnetoelectric coupling. Therefore, for simplicity we stick to the Pf theory. 

As in the magnetic case, the depolarizing energy (equal to ^q 1 P 2 for a thin ferroelectric 
film but less than this value for a heterostructure) is ignored in Eq. ([3]). It will be commented 
on later. Also, the gradient energy is ignored as the electric polarization is assumed to be 
aligned within each film. In reality, the polarization may be greatly changed near an interface 
but these more complicated effects are ignored here. 

The elastic energy density for both cubic phases is given for simplicity by an expansion 
in strain up to second order terms [22J: 



= \c^{u x f + u y f + uif) 
+\ctf (u x f + u y f + uif) 

+Ci$ + rigug + tiSJuW) ; (4) 

where p = m,p denotes the two phases and ufj are the strains. are components of the 
elastic compliance tensor for each material. 

The magnetostrictive energy density coupling the strains and magnetization in a cubic 



material is given by 
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FpM-eiO;,^) = B i (Ox) 2 ^ + (ay) 2 Uy7 ] + (a z ) 2 u^ 

I JD ( ( m ) I ( m ) I ( m )\ /r\ 

-r-L>2 [Q'xQ'yUxy T 0,yd z Uy Z T 0, z (l x Uzx J , ) 

where B l and B 2 , following the notation of Kittel, are the magnetoelastic constants with 



units Jm 3 and can be related to measured magnetostriction constants Ai o and A m 
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The energy density coupling strains and electric polarization is given similarly by [25] 

+^12 KAPy + P 2 Z ) + U yy {P 2 x + Pi) + U zz (Pl + P 2 y )) 

+q M {P x P y u%> + P y P z uf z + P a P a v£>) , (6) 

where P is the electric polarization and q n , q\ 2 and g 44 are the coupling coefficients with 
units JmC~ 2 . 

The total elastic contribution to the energy density of the effective medium can be found 
by adding together the correctly weighted energy densities, Eqs. (SJ)- (JH1) - 

Ptot.el. (Oj , Pj , U { ™ } , U W ) = — — I d m (.PpM-el , uj£° ) + F el ) ) 

+d p (F FE _ el (P u u^) + F el (u^))}. (7) 
The thicknesses of the magnetostrictive and piezoelectric films are given by d m and d p . 

B. Strain boundary conditions 

In order to find an approximate magnetoelectric coupling, we impose boundary con- 
ditions on the strains in the two respective phases and then minimize Eq. ([7]) in order to 
eliminate strain from the expression. The problem then is which boundary conditions to use. 
Harshe et al. considered four different cases of strain boundary conditions when deriving 
magnetoelectric coupling between magnetic and electric fields, H and E, in magnetostric- 
tive/piezoelectric thin film composites [loj]. Here only one set of boundary conditions will 
be considered. 

It is assumed that the heterostructure is mechanically clamped in the ^-direction, per- 
pendicular to the plane of each interface. This is illustrated in Fig. Ufa). Therefore the total 
distortion along the z direction must be zero and the boundary condition is jlo| 

d m u^ + d p u% = 0. (8) 

Furthermore, it is assumed that each interface has perfect mechanical bonding and therefore 
any distortion in the transverse directions will be equal in both phases: 

ui,^ = u^j. (9) 



"xx 
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FIG. 1: (Color online) Panel (a) shows a schematic of the heterostructure geometry. The structure 

piezoelectric (p) 



10( ] . Two repeats 



is clamped, with no friction, in the z-direction. The magnetostrictive (m) anc 
phases are mechanically bonded. This geometry corresponds to Case IV in Ref. 
only are drawn for simplicity. Infinite repeats of the magnetostrictive/piezoelectric thin films are 
considered in the effective medium calculation in the next section. In panel (b), a schematic 
drawing shows how the shear strains u yz may be different in different films. 



This ignores any mismatch strain between the two materials due to different unit cell sizes. 
This can be accounted for using methods used previously for piezoelectric films 



and 



would merely alter the boundary conditions by a constant while still ultimately giving rise 
to a magnetoelectric coupling. 

The shear strains correspond to a change in the angle between the unstrained cartesian 
coordinates. Therefore, the shear in the x — y plane must be equal in both films due to the 
perfect mechanical bonding: 

<? = «g. (11) 

However, the shears in the x — z and y — z planes have no restrictions. It is possible for 
a shear in the magnetostrictive material to be different than the shear in the piezoelectric 
material, as long as the height of the heterostructure in the z direction is conserved. This 
is illustrated schematically in Fig. [0(b) . 
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C. Solutions for the strains 



Substituting Eqs. ()8])-( TTTj) into Eq. ([7]), Uxx, Uyy , uiz and Uxy can be eliminated. Then 
solving 9 ^ n) = = x,y,z) and -^j = and -^y = 0, the strains that minimize the 
free energy are found. The simplest strain components to derive (since they are uncoupled 
from the other strains) are the shear components: 



_d m B2a x a y + d v q^P x P y 
B 2 a x a z 



c 44 



(m) 



(13) 



«2 = -"-^ (14) 



u[ m) 



L/ 44 

B 2 a y a z 



(15) 



(p) = _f44^ jT » QgN 



V z r<(m) 

lip) 
-"44 

By the symmetry of the boundary conditions, the transverse strains u xx and *Uyy cLI*6 
symmetric on exchange of subscripts x and y, whereas v,™ = —(d p /d m )uzJ is different. The 
transverse strains are therefore given by 

_ d m Bifi(ai) + d p qi 1 f 2 (P j ) + d p q 12 f 3 (Pj) . . 

C 1 C 2 

M OT = U xxiflx a yi Px ^ Py) (1^) 
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where functions /i, ji and f% are defined 



and 



—a 



(^mC'l ^ + ^pCll^X^P^ll ^ + ^mCi?) 



(d m C[ 2 ^ + dpCi2)(d p Cii 1 + <i m C^ 



v(p) 



(m) 



d m d p (C[ 2 ^ Ci2^) 
-o^ p (C (p) 



12 



(m)N 
12 , 



i- 

t(p) 



(m) 



+d m d p (C[ 2 ^ — C[ 2 ^) 2 



+P 



(d m C[ 2 + ^pCiIOl^p^ii ) ~^ d m C\i) 



y(p) 



(m) 



dmdpiCyz ^ C12O 



G>)\2 



12 



12 , 



1- 



/s(^) 



P 



^m(Cn^ + ^12 ^ — Cl2^)(Mm — ^p]Cl2 ^ + 2d p C 

+c[T\d m [d p - d m \c^ + [t4 + dl\c{ 2 )] 



12 ; 



Py Cq — P 2 (d m Cii + dpC^ ')ci 



y(m)N 



Ci = 



Co = 



~ a m L/ ll l U ll + U 12 ~~ U 12 J ~~ a m«pWl + «p U ll U ll 

~\~d m (— dpC'fi + [c[ 2 ^ — [— dpC'fi + (d m + d p )C[ 2 

"mV^ll — °12 I tt pl°ll ~~ °12 J' 

x (d m {c[i ^ + c| 2 ^) + ^(cjf + Ci2^) 



-2d m d p (c[ 2 ^ — c[ p ^ ' 



'12 



The zz-component is given by 
d 



u 



(m) 



= -T {- K + "»;) ^ f-Ci? + Ci?) - mlB.dpCs 
c 2 < v / 

+ [Pi + - C<?) + g 12 c 3 ) 

-P 2 2 (g 11 c 3 -2rf p g 12 (^ ) -^ ) ))} 
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where 



c 3 = d m (Ci? ] + Cg>) + d p (cS + C$). (26) 
D. Resulting free energy 

Substituting Eqs. f|T2l) - fl26l) for the strain back into the total strain free energy [Eq. ©], 
we obtain a energy density which takes the form 

Ftot.eiXdi, Pj) = X u af + X^a 2 a 2 + Y u Pf 

+Z ij a 2 i P* + Z: y a x P x a y P y , (27) 

where Einstein summation over variables i,j = x,y,z is assumed and % ^ j. The X coef- 
ficients (units Jm -3 ) represent the effect of strain on the magnetization, the Y coefficients 
(units JC _4 m 5 ) represent the effect of strain on the electric polarization, and the Z coef- 
ficients (units JC _2 m) give the strength of the magnetoelectric coupling. This represents 
many more terms than is calculated using simpler Landau expansions, where only the com- 
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4l4|. 



plete free energy of one of the two films is considered 

Due to the asymmetric boundary conditions, X xx = X yy ^ X zz , X xz = X yz ^ X xy , and 
so on for the other coefficients. The independent magnetic coefficients are given by 

B 2 d 2 c 2 

X xx = ~~Z7~-. ~tt 2 (^) 

2(u m -|- dp)c^C2 

d d B 2 r 
1(d m -\- dp)c2 

x d ™ (b 2 A \b\ \ 

XV (d m + d p )[c 1 c 2 dm C^+d p cS) [ ) 

A R2//i(p) _n( m )\ A R2 

(d m + d p )c 2 2(d m + d p )C[™ y 



where 



C4 — d m d p C[i + d m d p (c^± — (c[ 2 ^ ~ C*i2^) J 

Hd 2 m + d 2 p )CiT ] cff (32) 
r 2 - -d d (r {m) - r ,(p h 2 

-\-{d m C^ + dpC^ ^)(d m C^ + d p C^). (33) 
11 



The independent electric coefficients are given by: 



Y, 



y. 



-Qu d p4 + 2q 1 iq 12 d 2 p ci + qj 2 d p c 2 7 
dpjc^l 

dp (d m qfiC 3 — 4c? m c? p gn9i2(C^ — C{ 2 )) + ^d p qf 2 {d m C^i + d p C 

2 (d m + d v )c\ 



{my 
11 , 



Y, 



dpQn c 5 



die 



pHA4 
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(d m + d p ) Cl c| 2 (rf m + dp) (d m C { ™ ] + dpC^ 
d P ql 2 c\ , d 2 p q n qi2cl 



+ 



(cf m + d p )c l c 2 [d m + dp)cic; 



(<i m + dp)c 2 



| O^m ^p9ll(Cl2^ ~ Cl2 ^) — <^p<?12 



^m(ClV "I - 2C*12 ^ 2C^2 ) dpC 



~ <?n<?i2 [(eg ^ + eg + (2CiV + eg ^)d m dp 
dpq^ 

2C44 (<^m + <^p) 



where in Eq. (J34J) we have constants 



<4 ((eg 1 ) - c£°)(cg° - eg)) + effeg*) 
+«Mp (cg°(3cg> - cg>) + 2eg)(eg) - eg)) + eg^eg 1 ) - a 



12 , 



V™ (efr^eg) - seg)) + cfi^cg*' - 2c 



(m)> 
11 , 



+^(^ + 4)(e 1 ( 2 n)2 -eJr )2 )- 



"p !! ^11 



-2d m OT(e^ + er-e 



12 
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The magnetoelectric coefficients, which are of most interest, are given by: 

{dm 4~ dp)c^c^ 



d m dpB\ (^QnC 3 + 2q 12 dp{C{ 2 ^ ~ C12 



z ** = — n 1 3 (41) 



d m dpBi(q u (% + q 12 q 



Z xy = mp f Z ,\1 2 (42) 



d m d p B\ {^q\\d m {C\ 2 ■* — C^) 4- qi 2 {d m C[^ 4" rfpCn ' 

^ xz = 7~j i — , x 2 = ^j/z (43) 

(a m + dp)c 2 

d m d p B l [qndpiC^ - C$) + q 12 c 3 ) 

Z zx = ^ — - 2 '- = Z zy (44) 

[dm 4" d p )c 2 



Z* = ^44^ (45) 



d m d v q^B 2 
{d m 4" dp) {d m C 44 -|- dpC ^ 

Notice that if the thickness of either films vanishes (d m / p — >■ 0), then the magnetoelectric 
coupling vanishes. There is also an optimum ratio d m /d p such that the magnetoelectric 
coupling is maximized, as will be shown in Section IHII 

Now that strain has been eliminated, the effective free energy due to strain [i*tot.ei> 
Eq. (|27|) ] can be added to the unstrained free energies of the piezoelectric and magnetostric- 
tive phases [Ffm 4- Ffe, Eqs. (j2j) and (jHj)]. Before calculating the ferromagnetic resonant 
frequency of such a system in Section IHII we first estimate how the strain alters the mag- 
netic and piezoelectric constants {K and a) and estimate the strength of the magnetoelectric 
coupling for a BaTi0 3 /NiFe 2 4 heterostructure. 



E. BaTi03/NiFe204 coupling strength 

A heterostructure comprising equal thickness films {d m = d p ) of piezoelectric BaTiC>3 and 
magnetostrictive NiFeoC^ is considered. The ferroelectric parameters for BaTiC>3 in Eq. ([3]) 
are taken from Ref. [20| and are: a x = -2.772 x 10 7 J rn (T 2 , a u = -6.476 x 10 8 J m 5 CT 4 , 
a 12 = 3.230 x 10 8 J m 5 C" 4 , a m = 8.004 x 10 9 J m 9 C" 6 , a ll2 = 4.470 x 10 9 J m 9 C~ 6 
and 0123 = 4.910 x 10 9 J m 9 C -6 . The piezoelectric and elastic compliance parameters 
are also taken from Ref. [20]: q u = -14.20 x 10 9 J m C~ 2 , q l2 = 0.74 x 10 9 J m C~ 2 , 
g 44 = -3.14 x 10 9 J m C" 2 , Cf? = 27.50 x 10 10 J m" 3 , C$ = 17.90 x 10 10 J m" 3 and 
C$ = 5.43 x 10 10 J m" 3 . 
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The parameters for NiFe 2 4 are taken from Ref. Q: = 21.99 x 10 10 N m" 2 , C[™ ] = 
10.94 x 10 10 N m" 2 , C { ™ ] = 8.12 x 10 10 N m" 2 and /i M = 0.32 T (4ttM = 3200 G 
in CGS). The demagnetizing ener gy density is then 8.2 x 10 4 J m -3 . For the anisotropy 
constant the value found by Weisz 26j is used, namely K = 4 x 10 3 J m -3 . An applied field 
HoH = 0.1 T x is also applied to ensure that a preferred direction of magnetization exists 
and later the equations of motion will be well-defined. 

With these material parameters, the equilibrium magnetization and electric polarization 
can be found by calculating the free energy minima for each possible direction of the order 
parameters [181 ] . The equilibrium corresponds to the state with the lowest energy minima. 
The spontaneous electric polarization and magnetization are found to both lie along the 
x axis (see Fig. CD), ie. 

a = (l,a y ,a z ) (46) 
P = (P + p x ,p v ,p g ), (47) 

where the lower case symbols a y , a z « 1 and Pi << P (i = x, y, z) represent small dynamic 
quantities that will be relevant in Section [TTT1 This result is obtained whether we consider the 
demagnetizing and depolarizing energy to be that of a thin film or that of bulk. In reality, 
the demagnetizing and depolarizing energy is somewhere between these two limits for the 
heterostructure and therefore the result holds for this geometry also. Then, by ignoring at 
first the contribution to the free energy due to strain [Eq. ([27}) ] and only considering Eq. (j3J), 



the polarization in zero applied field is 17] 



P = 11 v 11 — 1 111 = 0.265 Cm- 2 . (48) 
3a m 

Now it is possible to estimate the effect that the strains have on the magnetic and electric 
systems. Firstly, the strains alter an by over 10%: a\ x — > (an + Y xx ) and a\ x — > (an + Y zz ), 
where Y xx = -1.66 x 10 8 J C' 4 m 5 and Y zz = -1.03 x 10 8 J C" 4 m 5 , using Eqs. flU} and 
( 135]) . This in turn alters the value of the polarization calculated without strains in Eq. ( 148]) 
by 8% (P = 0.286 C m- 2 ). 

All of the strain-mediated magnetic, electric and magnetoelectric coefficients are listed 
in Table |I] for the material parameters listed above. The electric coefficients are multiplied 
by P 4 = 0.286 4 and the magnetoelectric coefficients are multiplied by P 2 = 0.286 2 so that 
every term has units of energy density and their strengths can be compared. 
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TABLE I: Values of the strain- mediated magnetic (Xij), electric (Yij) and magnetoelectric (Zij) 
coefficients for a BaTi03/NiFe2 04 thin film heterostructure with d m = d p , derived from Eqs. ()28|) - 
(|45|) . coefficients are multiplied by P 4 = 0.286 4 C 4 m~ 8 and coefficients are multiplied by 
P 2 so that all values are in units J m~ 3 . 



X xx 




X X y 


X XZ 


-9.86 


-6.64 


-39.7 


-86.4 


YxxP 


Y ZZ P 4 


YxyP 4 


Y XZ P 


-1.11 x 10 6 


-7.14 x 10 5 


1.27 x 10 6 


-6.88 xlO 4 


ZxxP 


Z ZZ P 2 


ZxyP 


Z X zP 2 


6.61 x 10 3 


-4.35 x 10 3 


-3.99 x 10 3 


-524 


z zx p 2 


^xy r 






587 


4.99 x 10 3 







Strain has a large effect on the ferroelectric parameters, as illustrated by the 8% change in 
P, but a very small effect on the magnetic parameters. This is perhaps not surprising since 
lattice displacements (strains) are intimately related to ferroelectricity but are not necessary 
for magnetic ordering. For example, the magnetic coefficients Xa are all on the order of 1 - 
80 J m -3 (see Table H]) and so strain changes the anisotropy energy by under 1%. From now 
on we therefore ignore the effect of the Xa . 

The fact that strain has a large effect on the piezoelectric phase and a relatively small 
effect on the magnetostrictive phase leads to a strain-mediated magnetoelectric coupling 
which has an intermediate strength. It is proportional to a product of magnetostric- 
tive/piezoelectric coefficients, such as B\qu. For example, Z XX P 2 = 6.61 x 10 3 J m~ 3 . 
While this value represents less than 0.3% of the ferroelectric energy density contributions 
(characterized by a±), it is larger than the magnetic anisotropy energy in an unstrained 
system. Therefore, we can expect in our subsequent calculation for the high frequency sus- 
ceptibility that the ferromagnetic resonance (FMR) frequency will be shifted by an electric 
field while the affect of a magnetic field on the electric modes will be negligible. 
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III. HIGH FREQUENCY EFFECTS OF MAGNETOELECTRIC COUPLING 



To calculate the high frequency susceptibility x °f the heterostructure, the effective 



medium method presented in Ref. |3| is used. Unlike in Ref. [3j, the full 6x6 suscepti- 
bility is not solved for. Instead, only xTz = m fV h e z s is solved for to demonstrate the action 
of the magnetoelectric coupling on the FMR frequency. In Section [Til the number of repeats 
of the thin films is considered finite, as illustrated in Fig. [TJ However, the effective medium 
method considers the number of repeats to be infinite. The results for the frequency should 
be accurate for a large number of unit cell repeats and will be qualitatively correct for a 
small number of repeats. 



A. Equations for the resonant frequencies 

The calculation involves solving the linearized magnetization and electric polarization 
equations of motion together with Maxwell's boundary conditions for the dipolar fields. 
The total energy density F = Ffm + -Ffe + ^tot.ei., the sum of Eqs. © an d (127)1 . is 
substituted into the equations of motion (in SI units): 

da ( 1 5F\ 

-JI = 70 * -T7TZ (49) 



dt V M 5a 

d 2 P JF 



it* ~ e ° f 6P- (50) 

where 7 is the gyromagnetic ratio (7 = 27r2.89 x 10 10 rad Hz T _1 ) and / is the effective 
mass term (/ = (2tt) 2 1.5 x 10 26 Hz 2 ). Then the equations are linearized using Eqs. (|46p and 
(|47j) and small oscillatory solutions are assumed (a y / z oc e~ l0Jt and Pi oc e~ luJt ). The resulting 
equations of motion for the magnetization in the magnetostrictive phase are: 

- l -^a y = = a z [B K + B - B XX (P) - B XZ (P)] 

hW (51) 



iuj 



^ = -^ + B -B U P) + B xy(P)] 

Z:yP Py + ^h v m \ (52) 



M 



where B = [IqH represents the static magnetic field along the +x direction, B K = is 
the effective anisotropy field, Bij(P) = ^jf- is the effective magnetic field due to magne- 
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toelectric coupling and hz and h™ are components of the dipolar magnetic field in the 
magnetostrictive material. 

In Eq. (152]) . it can be seen that the magnetoelectric coupling energy Z* y a x a y P x P y leads to 
a dynamic coupling between the magnetization and the ^-component of the electric polar- 
ization in the piezoelectric phase. Moreover, the other magnetoelectric coupling terms alter 
the effective static magnetic field felt by the system and therefore will alter the resonant 
frequency. 

Since p y is dynamically coupled to the magnetization equations of motion, we also write 
its equation here: 

U) 2 

7Py = ~2Py [&1 + ("12 + ^y)-P 2 + «112-P 4 ] ~ Z* Pa y + S> ( 53 ) 

where e y is a component of the electric dipolar field. 

Maxwell's boundary conditions between the films for the magnetic dipolar fields in the 
magnetostatic limit are that the in-plane component of the magnetic field and the out-of- 
plane component of the magnetic induction must be continuous, ie. 

/4 m) + a z M = h£> = — (54) 

*5S = C < 55 > 

where C is a constant, defined for convenience in the calculation and h^ p ' is the magnetic 
dipolar field in the piezoelectric films. The susceptibility components due to the out-of-plane 
dipolar field jif 1 ^^ must be calculated first to properly include the dipolar field effects due 
to the boundary conditions. It is this step of including Maxwell's boundary conditions before 
solving for the equations of motion that sets this type of effective medium method apart 
from other effective medium methods and which leads to a more sophisticated model for the 
demagnetizing effects on the frequencies. 

The zz component of the magnetic susceptibility of the effective medium can be calculated 
by setting h y = = e y and by using the weighted average: 

v m = ™? = d m a z M 

Xzz ~hf dm h^ + d p h^- [ j 

Eqs. (15T|) - (154|) can be rearranged so that a z , a y , p y , h!^ 1 and h z v ^ are all written in terms of 
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the constant C, which therefore cancels in Eq. (|56l) . The result is 

m dmHoMn? (go + B K - B xx + B xy + jgggy 7 
Xzz d m {u*-ul) + d v {uj*-uj) 

where Ub is the bulk ferromagnetic frequency under strain given by: 

^ ( B +B B +B 4- fe ° Z °y P2 



7 2 V ^ 2 - <^ 

x (B + B K - B xx - B xz + fi M) , (58) 

Uf is the thin film ferromagnetic frequency under strain given by: 

/ ftoZfyP 2 

— = \B + B K — B xx + B xy + 



x (B + fi x - 5 ra - B xz ) , ' (59) 

and ut is the frequency of the transverse phonon mode associated with the ferroelectric 
polarization, with oscillation in the y direction 

c4 = fe Q 2 (at + (a 12 + Y xy )P 2 + a U2 P 4 ) . (60) 



It can be noticed that Eqs. (158]) and p9j) , are implicit equations in to. However, it can 
be shown that ojb/f (typically in the low GHz) and ujt (typically in the low THz regime) are 
very widely spaced in frequency. Therefore, the magnetic frequencies can be approximated 
very well by ignoring the term proportional to Z* y in Eqs. fl58l) and (159]) . 

The pole of Eq. ( 157]) is given by the positive solution to 

= d m (u 2 - u 2 ) + d p (u 2 - u 2 ) (61) 

and corresponds to the ferromagnetic resonance frequency. This equation shows how the 
effective medium method calculates a resonant frequency for the heterostructure which is 
between that of bulk Ub and that of a thin film ojf. In Ref. 3J it was shown that the effective 
medium FMR frequency reduces to the correct results in these well-known limits. 

B. Electric field shift of the FMR frequency 

Since applying an electric field alters the effective magnetic fields B XX (P) and B XZ (P) 
in Eqs. (158]) and (159]) . it shifts the FMR frequency of the effective medium. This is shown 
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in Fig. [2j In panel (a) xTz (Eq- (156]) ) is plotted as a function of frequency for no applied 
electric field (solid line) and for an applied electric field along the +x direction equal to 
10 7 V/m (dashed line). The pole in the susceptibility corresponds to the FMR frequency. 
Since damping is not included in the calculation, the peaks are not rounded at all, as is seen 
in real experiments. The FMR frequency is plotted as a function of electric field in panel 
(b). 
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FIG. 2: (Color online) The zz-component of the magnetic susceptibility as a function of frequency 



(panel (a)) for a BaTi03/NiFe20"4 heterostructure, with each film of equal thickness d r , 



d n . A 



peak occurs at the ferromagnetic resonance (FMR) frequency. The peak is shifted downwards by 
application of an electric field in the positive x-direction (dashed line). In panel (b) the FMR 
frequency is plotted as a function of applied electric field strength. 



Although it is encouraging to see a shift in the FMR frequency when an electric field is 
applied, a large field is needed to see a 0.2 GHz shift and this appears too small for signal 
processing applications. However, it is slightly larger than FMR shifts seen in most recent 



experiments on similar systems, for example, Refs. 27|, |28[ . Alternatively, experiments may 
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apply a constant driving frequency and sweep the applied magnetic field H and then the 
resonant field shift as a function of E may be measured. At a constant driving filed of 9 GHz, 
our calculation gives a shift in H of 350 Oe (hqAH = 0.35 T) when a field E = 10 s V/m 
is applied compared with when there is no field. This value is higher than measurements 
made, for example, by Bichurin et al. [2J on ferrite/PZT samples but on the same order of 
magnitude as that measued by Liu et al. 29]. However, even larger effects have recently 



been seen when the ferrite material is replaced with the amorphous alloy FeGaB 
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FIG. 3: (Color online) The zz-component of the magnetic susceptibility as a function of frequency 
for a BaTi03/NiFe204 heterostructure (d m = d p ) at high frequencies. A peak occurs at the 
transverse phonon resonance. The peak is shifted upwards by application of a positive electric field 
(dashed line). Notice that the vertical scale is 10 4 times smaller than that of Fig. [2ja). 



It should be noted that there is a second, weak pole in the plot of xTz as a function of 
frequency outside the range shown in Fig. |2£a). This occurs at the transverse phonon mode 
frequency and is due to the dynamic magnetoelectric coupling discussed. It could not be 
predicted in an EMM where the polarization dynamics were ignored. Since it is an electric 
mode frequency, it is highly tunable with electric fields, as illustrated in Fig. |3j A field of 
E = 10 7 V/m shifts the frequency upward by roughly 32 GHz. Although it is highly tunable, 
this resonance is too weak and too high to be useful for microwave applications. 

To try to match with particular experiments, the exact material parameters and geom- 
etry and strain boundary conditions must be input into the model. In particular, most 
experiments are performed with the poling direction out-of-plane whereas the results in 
this section have the polarization and applied electric field in-plane to simplify the effective 
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medium method. Changing the poling direction to out-of-plane (therefore able to apply 
larger electric fields) and modeling specific experiments represents an exciting extension to 
the work presented here. 

Only one other set of strain boundary conditions is considered to see what effect different 
boundary conditions have on the solutions demonstrated in Figs. [2] and |3J The mechanical 
clamping in the ^-direction [Eq. flH])] is relaxed and Ui£ and 11?} are decoupled. Repeating 
the entire calculation for the strain-mediated magnetoelectric coupling and the resulting 
high frequency dynamics, it is found that the change in boundary conditions has little effect 
on the FMR frequency, changing the values calculated by under 0.01 GHz. However, the 
frequency of the transverse phonon mode, ut, is shifted down in frequency by roughly 
10 GHz. This reflects again the fact that the piezoelectric system is sensitive to changes in 
the strain whereas the magnetic system is weakly affected. Also, for this geometry, Z xx is 
the magnetoelectric coupling term which most affects the FMR frequency (see Eqs. (|58|) and 
( |59l ) and this term is relatively unchanged when the u zz strains are altered. Changing the 
boundary conditions for in-plane strains may therefore alter the FMR frequency by a larger 
amount. 

C. Thickness effects 

Before concluding this section, the effect of the relative film thicknesses will be dis- 
cussed. The ratio of the thickness of the piezoelectric versus the magnetostrictive film 
(d p :d m ) changes both the magnetoelectric coupling strength and the magnetic dipolar fields 
in the heterostructure. Therefore the effective medium calculation for the FMR frequency is 
altered through Eq. floTj) . The result of these two effects is shown in Fig. HI where the FMR 
frequency is plotted as a function of the volume fraction of magnetostrictive material. The 
solid line shows the result in zero applied electric field and the dashed line shows the result 
when E = 10 7 V/m. The dotted line at the top shows the result of the effective medium 
calculation if the effect of strains is ignored (By = 0). 

It can be seen in Fig. H] that the strain coupling between the magnetostrictive and 
piezoelectric films changes the FMR frequency significantly. For a volume fraction of 
d m / (d m + dp) = 0.4, the FMR frequency in the strained system (solid line) is half that 
for the unstrained system (dotted line). Fig. H] also reveals that a volume fraction of around 
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FIG. 4: (Color online) The FMR frequency, calculated using Eq. (|61h . as a function of the magne- 
tostrictive material's volume fraction d m /{d m + d p ). The solid line shows the result in zero applied 
electric field and the dashed line shows the result when E = 10 V/m. The dotted line shows the 
result if there were no strain coupling to the piezoelectric. 

0.5 results in the largest possible change in the FMR frequency on application of an ap- 
plied electric field, although this effect is weak. All three frequency calculations agree in 
the limits d rn — >■ and d p — > since then the strain coupling (and electric field) have no 
influence on the FMR frequency. For d m — > (volume fraction is 0), the frequency of a thin 
magnetostrictive film is found. For d p — > (volume fraction is 1), the frequency of the bulk 
magnetostrictive material is found. In between, it can be seen how the effective medium 
method finds frequencies between these two limits since the dipolar energy is between these 
two extremes. 

IV. CONCLUSION 

In conclusion, it has been shown how strain-mediated magnetoelectric coupling energies 
can be determined for magnetostrictive/piezoelectric thin film heterostructures. This is done 
using Landau free energy expansions in terms of the magnetization, polarization and strains 
of both thin films, and then eliminating strain using boundary conditions. 

The resulting terms are used to calculate the electric field shift in the FMR frequency 
in a BaTi03/NiFe204 heterostructure. A shift of 0.2 GHz is calculated for an applied 
field E = 10 7 V/m. This is on the same order of magnitude as some experiments. The 
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particular geometry and strain boundary conditions that are assumed may affect the strength 
of the derived magnetoelectric coupling terms and therefore affect the shift in the FMR 
frequency. An analysis of different materials will enable predictions to be made on which 
heterostructures are best for possible applications. 

What has not been mentioned so far in this work is the fact that magnetoelectric (fer- 
rimagnetic/ferroelectric) heterostructures in which there is little strain-mediated coupling 
between magnetic and electric phases have also been shown to have a large magnetoelectric 
coupling and shifts in the FMR frequency may occur of the same magnitude as in mechan- 



ically bonded systems [3J 



331 ] . This is due to the intrinsic hybridization of electric and 



magnetic fields and is not captured in our theory since the magnetostatic and electrostatic 
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34j may be used to calculate 



limit is assumed. Methods such as those described in Refs. 
the FMR frequency beyond the magnetostatic limit. Future work will involve adding strain- 
mediated magnetoelectric couplings, like those derived in this paper, into such a theory to 
incorporate both the effects of strains and hybridized fields and to determine their relative 
importance. 
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